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We study the QCD phase diagram, in particular we study the critical points of the two main 
QCD phase transitions, confinement and chiral symmetry breaking. Confinement drives chiral 
symmetry breaking, and, due to the finite quark mass, at small density both transitions are 
a crossover, while they are a first or second order phase transition in large density. We study 
the QCD phase diagram with a quark potential model including both confinement and chiral 
symmetry. This formalism, in the Coulomb gauge hamiltonian formalism of QCD, is presently 
the only one able to microscopically include both a quark-antiquark confining potential and 
a vacuum condensate of quark-antiquark pairs. This model is able to address all the excited 
hadrons, and chiral symmetry breaking, at the same token. Our order parameters are the 
Polyakov loop and the quark mass gap. The confining potential is extracted from the Lattice 
QCD data of the Bielefeld group. We address how the quark masses affect the critical point 
location in the phase diagram. 



1 Introduction 

Our main motivation is to contribute to understand the QCD phase diagram for finite T and 
H, to be studied at LHC, RHIC and FAIR-'-. Moreover, our formahsm, in the Coulomb gauge 
hamiltonian formalism of QCD, is presently the only one able to microscopically include both 
a quark-antiquark confining potential and a vacuum condensate of quark-antiquark pairs. Thus 
the present work, not only addresses the QCD phase diagram, but it also constitutes the first 
step to allow us in the future to compute the spectrum of any hadron to finite T . 

Here we address the finite temperature string tension, the quark mass gap for a finite current 
quark mass and temperature, and the deconfinement and chiral restoration crossovers. We 
conclude on the separation of the critical point for for chiral symmetry restoration from the 
critical point for deconfinement. 



2 Fits for the finite T string tension from the Lattice QCD free energy 

At vanishing temperature T = 0, the confinement, modelled by a string, is dominant at moderate 
distances, 

V{r) c^^ + Vo + ar . (1) 

At short distances we have the Luscher or Nambu-Gotto Coulomb due to the string vibration 
+ the OGE coulomb, however the Coulomb is not important for chiral symmetry breaking. At 
finite temperature the string tension (t{T) should also dominate chiral symmetry breaking, and 
thus one of our crucial steps here is the fit of the string tension a{T) from the Lattice QCD data 
of the Bielefeld Lattice QCD group, ISEIlEll. 




The Polyakov loop is a gluonic path, closed in the imaginary time (proportional to the 
inverse temperature T~^) direction in QCD discretized in a periodic boundary euclidian Lattice. 
It measures the free energy F of one or more static quarks, 

P(0) = iVe-^"/^ , P"(0)P^(r) = Are-^^?(^)/^ . (2) 

If we consider a single solitary quark in the universe, in the confining phase, his string will travel 
as far as needed to connect the quark to an antiquark, resulting in an infinite energy F. Thus 
the 1 quark Polyakov loop P is a frequently used order parameter for deconfinement. With the 
string tension cr(T) extracted from the qq pair of Polyakov loops we can also estimate the 1 
quark Polyakov loop P{0). At finite T, we use as thermodynamic potentials the free energy Fi 
and the internal energy Ui , computed in Lattice QCD with the Polyakov looplSEHEE] 

They are 

related to the static potential V{r) = —fdr with Fl{r) = —fdr — SdT adequate for isothermic 
transformations. In Fig. [T] we extract the string tensions ct(T) from the free energy Fi(T) 
computed by the Bielefeld group, and we also include string tensions previously computed by 
the Bielefeld group^. 

We also find an ansatz for the string tension curve, among the order parameter curves of 
other physical systems related to confinement, i. e. in ferromagnetic materials, in the Ising 
model, in superconductors either in the BCS model or in the Ginzburg-Landau model, or in 
string models, to suggest ansatze for the string tension curve. We find that the order parameter 
curve that best fits our string tension curve is the spontaneous magnetization of a ferromagnet 
^, solution of the algebraic equation, 

tanh . (3) 
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In Fig. [1] we show the solution of Eq. [3] obtained with the fixed point expansion, and compare 
it with the string tensions computed from lattice QCD data. 

3 The mass gap equation with finite temperature and finite current quark mass. 

Now, the critical point occurs when the phase transition changes to a crossover, and the crossover 
in QCD is produced by the finite current quark mass mO, since it affects the order parameters 
P or fj, and the mass gap m(0) or the quark condensate {qq). The mass gap equation at the 
ladder /rainbow truncation of Coulomb Gauge QCD in equal time reads, 

a dk lA{p,k,fj.)m{k)p - lB{p,k,fi)m{p)k 

-(^) = -o+^y^ ^ v^T^t^F ' ^ ^ 




Figure 2: (left) The mass gap m(0) solution of as a function of the quark current mass mo, in units of a = 1. 
(right) Sketch of the effect of mo on the crossover versus phase transition of choral restoration at finite T. 




Figure 3: Sketches of the saturation of confinement (left), and of the corresponding crossover in the order param- 
eter P polyakov loop (righ). 
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The mass gap equation dH) for the running mass m{p) is a non-hnear integral equation with 
a nasty cancellation of Infrared divergences ^ We devise a new method with a rational 

ansatz, and with relaxation, to get a maximum precision in the IR where the equation is nearly 
almost unstable. The solution m{p) is shown in Fig. [2] for a vanishing momentum p = 0. 

At finite T, one only has to change the string tension to the finite T string tension a{T)^^^, 
and also to replace an integral in w by a sum in Matsubara Frequencies. Both are equivalent to 
a reduction in the string tension, a ^ a* and thus all we have to do is to solve the mass gap 
equation in units of a* . The results are depicted in Fig. [2j Thus at vanishing niQ we have a 
chiral symmetry phase transition, and at finite mg we have a crossover, that gets weaker and 
weaker when mo increases. This is also sketched in Fig. [2j 
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4 Chiral symmetry and confinement crossovers with a finite current quark mass 

In what concerns confinement, the linear confining quark- antiquark potential saturates when the 
string breaks at the threshold for the creation of a quark-antiquark pair. Thus the free energy 
F(0) of a single static quark is not infinite, but is the energy of the string saturation, of the 



order of the mass of a meson i. e. of 2mo. For the Polyakov loop we get, 



P(0) ~ iVe-2™o/T ^ 

Thus at infinite mo we have a confining phase transition, while at finite rriQ we have a crossover, 
that gets weaker and weaker when niQ decreases. This is sketched in Fig. [3l 

Since the finite current quark mass affects in opposite ways the crossover for confinement 
and the one for chiral symmetry, we conjecture that at finite T and /u there are not only one but 
two critical points (a point where a crossover separates from a phase transition). Since for the 
light u and d quarks the current mass mo is small, we expect the crossover for chiral symmetry 
restoration critical to be closer to the ^ = vertical axis, and the crossover for deconfinement 
to go deeper into the finite fi region of the typical critical curve in the QCD phase diagram 
presently assume 
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